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ABSTRACT 


A gas-cooled, cylindrical, nuclear reactor is used as the 





basis in the theoretical derivation of the coolant's enthalpy increase 
and pressure loss acrose the reactor. Turbulent flow is assumed in 
the coolant passages of the reactor, and the Reynolds analogy is 
used as the basis for correlating the heat transfer coefficient and 
the friction coefficient. 

The general equations are cerived and two examples, a 
nuclear-hydrogen rocket and a ramjet, are given to demonstrate 


applications of the general results. 





SUMMARY 


The development of nuclear reactors has led to a new class of heat 
transfer problems. Heat is produced within a solid fuel at a rate that 
is almost independent of the coolant flow rate, at least for gas-cooled 
reactors. The power density varies with position in general, although 
it can be made nearly uniform with the introduction of appropriate re- 
flectors and with variation of fuel concentration. However, these 
measures are not desirable, particularly if there are weight restric- 
tions for the reflector and cost restrictions on fuel elements. 

Aside from considerations of weight and cost, the most important 
characteristic of a reactor is the ratio of pumping work to the increase 
of thermodynamic availability of the coolant passing through the re- 
actor. For any single coolant passage through the reactor, this ratio 
can be determined readily from the Reynolds analogy, for any speci- 
fied power distribution along the length of the channel, and a given 
coolant flow rate through the channel or a given pressure drop over the 
reactor. The limitation on the solid fuel temperature will fix the 
roinimum pressure drop or flow rate. If the power density is the same 
for all channels, the results are simply summed over the reactor 
cross-section. 

When the power density varies from channel to channel, there is 
no simple way of calculating pressure drop or temperature rise 
through the reactor. The pressure drop for all channels must be the 
same, but the flow rate and temperature rise will vary from channel to 
channel. Apparently the calculations of temperature and flow rate dis- 
tributions over the reactor have been done by numerical computation 


in the past. An iterative procedure is developed in this thesis for 





determination of the overall reactor performance in analytical form. 
The sole restriction on the procedure is that the Mach number of the 
flow entering any coolant channel be appreciably less than unity. This 
restriction is not a serious one because the pressure drop resulting 
from higher entering Mach numbers would be generally unacceptable. 

The principal objective of this thesis is the development of the 
method of calculating reactor performance for any specified power den- 
sity distribution. In final form the method gives the pressure drop, 
the temperature distribution and mass flow rate distribution in the 
coolant leaving the reactor, as well as maximum wall temperatures, 
where porosity, coolant channel geometry, and power density distribu- 
tion are prescribed. The temperature of the coolant leaving the re- 
actor can be made uniform over the reactor cross-section by introducing 
appropriate restrictions in the individual coolant channels. This pro- 
cedure increases the average temperature rise through the reactor for 
a given maximum wall temperature, but does so at the cost of in- 
creased pressure drop. Whether restriction of the coolant flow channels 
in this way is advantageous or not will depend upon the particular ap- 
plication. The method of analysis developed here can be used to deter- 
mine the optimum restriction of the flow channels for any specific 
example. 

To illustrate the application of the method for calculating reactor 
performance, the nuclear rocket and nuclear ramjet were chosen as 
examples. The specific impulse for the rocket is found to be higher 
when the coolant passages are restricted to give a uniform coolant out- 
let temperature than for no restriction. The thrust coefficient of the 


ramjet is not improved appreciably by restricting the flow rate. 
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SYMBOLS 


Coolant passage cross-sectional area 
Radius of reactor 

Friction coefficient 

Heat transfer coefficient (Stanton number) 
Specific heat at constant pressure 

Axial power distribution function 

Radial power distribution function 
Conductivity 

Coolant passage length 

Total mass flow rate through reactor 
Mass flow rate per unit area 

Coolant passage wetted perimeter 

Gas pressure (static) 

Power density i: reactor; heat produced per unit volume 
Maximum power density in reactor 

Gas constant 


Gas temperature (static, averaged over coolant passage 
cross-section) 


Gas velocity (averaged over coolant passage cross- 
section) 


Porosity; void volume/total volurne 

By? 2 
C,-~- , reat transfer parameter 

ie. 

C oe @ 
~“{£{ PL , friction parameter 
oa 
r/a@ , dimensionless radius from centerline of reactor 
Viscosity 


Kinematic viscosity 


x/£L , dimensionless distance frora coolant passage 
entrance 











Subscripts: 


oe 


SYMBOLS 


Gas density (averaged over coolant passage cross- 
section) 


Frandt! number 


Entrance to coolant passage 
Exit from coolant passage 
Wall of coolant passage 


Total 








I INTRODUCTION 


The power distribution in a nuclear reactor is proportional to 
the neutron density which, in a bare cylindrical reactor, is highest in 
the center and decreases both axially and radially away from the cen- 
ter. When the reactor is used as a source of heat, the variation in 
power density distribution results in an uneven heat supply to the 
coolant. 

The neutron density can be made more uniform by enclosing the 
reactor with a reflector which causes some of the leaking neutrons to 
diffuse back into the core, or by varying the fucl loading so that the 
fuel concentration is higher towards the edges than in the center. 
However, both of these methods have disadvantages. The reflector 
material is usually quite bulky and relatively heavy. As a result, in 
several applications of nuclear reactors as a source of power, spe- 
cifically, for airborne propulsive applications, it is necessary to 
restrict the amount of reflector material that might normally be used. 
Fuel concentration variation is undesirable since reactors are usually 
constructed by assembling many small lamina, or rods, to obtain the 
final configuration. For manufacturing ease and interchangeability, 
it is more desirable to have uniformity of the fuel elements. 

The uneven power density distribution of the nuclear reactor 
results in the coolant having temperature gradients in the exhaust if 
the mass flow in each of the coolant channels ia the same. The cool- 
ant in the channel through the center of the core will be heated the 
most, whereas the flow in the outer channels will be much cooler due 
to the lower heat generation at the reactor extremities. Flow of this 
type is usually undesirable, since the surface temperature of the chan- 


nel walls in the outer channels is much lower than the maximum 


ak 


allowable, and thus the overall efficiency of the power plant is reduced. 
In addition, the temperature gradients of the exhaust of this type flow 
lead to mixing losses in the nozzle. 

Theoretically, these exhaust gas temperature gradients can be 
eliminated by restricting the mass flow in the outer channels so that 
the mass flow in every channel is proportional to the heat release in 
that channel. The result is that the outer channels, where the heat 
release is less, are cooled to a lesser degree so that each channel 
operates at its maximum allowable wall suriacc temperature, and the 
temperature of the exhaust gases from the reactor is uniform. This 
mass flow restriction can be accomplished by varying the flow channel 
cross-sectional areas continuously across the reactor core so that for 
the same overall pressure drop across the channel the mags flow is 
smaller. The same effect can be obtained by installing orifices in the 
channels to induce a pressure drop so that again the mass flow is pro- 
portional to the heat release of the particular channel. Of these two 
methods, the former is preferable since it utilizes all the available 
pressure drop to improve the heat transfer coefficient. However, it 
would not usually be practical to assemble a reactor in this manner, 
and so the orificing technique is considered to be more practical. 
Elimination of the exit temperature gradients and increasing the wall 
surface temperatures of all channels to the maximum allowable by 
either of these methods increases the overall pressure drop across 
the reactor and reduces the overall mass flow, in addition to increas - 
ing the average temperature rise per unit mase flow. Because of the 
increased pressure losses and the reduction in mass flow, the over- 


all performance of the reactor must be evaluated to determine the 








performance gain, if any, by orificing tc ascertain if the increased 
design and assembly complications are warranted. 

W. B. Hall (1) has worked out an example of the restricted flow 
reactor configuration. The basic reactcr assumptions made by Hall 
and those of this thesis are similar. However, the analysis here is 
more general, including allowance for higher pressure drop and 
higher velocities in the coolant channels, and a procedure for finding 
overall performance when the coolant outlet temperature varies 

radially. These improvements in the analysis are needed for appli- 
cation to reactors for propulsion. 

In order to keep the general discussion of heat transfer and 
pressure drop in gas-cooled reactors for propulsive applications with- 
in reasonable limits, simplifying assumptions are required. For 
this report, the reactor was taken as cylindrical, with radial sym- 
metry in fuel loading and coolant passage distribution. The coolant 
passages were assumed to be of uniform size, and the cross-sectional 
area was assumed constant over the length of the reactor. The power 
density of the reactor was assumed to vary axially and radially ac- 
cording to simple algebraic laws. Since the analysis was made pri- 
marily to investigate the result of using a reactor as a type of heat 
exchanger, or source of heat, reactor physics as such waa not con- 
sidered, and no account was taken of the power distribution shift 
during operation, etc. The porosity of the reactor was not assumed 
constant in the derivation of the general cquations, but was assumed 
constant in specific examples worked out to demonstrate the use of 


the results. 
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The evaluation of the overall performance of a propulsion 
systern involves determining the overall enthalpy rise of the working 
fluid and the losses of the system. For a reactor of given design, 
the controlling factor in onerating power {is normally some maximum 
allowable temperature in the system. This limit may be set by the 
phase stability of the fuel elements, or the moderator, by the allow- 
able thermal stresses in the fuel or some other part of the system, by 
the influence of temperature on corrosion, or some other particular 
thermal effect. For a given design and this limiting temperature, 
the temperature distribution in the reactor core can be determined, 
and from this, the maximum allowable coolant wall temperature 
supiect to the limiting temperature of the particular configuration. 
Since this report was intended to be quite general, the controlling 
factor for reactor operation was selected as the maximum allowable 
coolant wall temperature which can be determined when given a par- 
ticular configuration with type of fuel, concentration, etc. Thus the 
problem to be solved involves calculating the overall enthalpy rise 
and pressure losses for a given maximum wall temperature and given 
initial conditions. 

The high power requirements of today's rockets and ramjets 
and the expanding technology in the nuclear reactor field make the pos- 
aibility of using the gas-cooled reactor as a source of propulsive 
power for these vehicics a distinct reality. Applications for rocket 
engines hold promise, since the effective exhaust velocity is propor- 
tional to the inverse of the square root of the molecular weight, which 


suggests the possibility of using low molecular weight gases such as 











~5 
hydrogen for the working fluid. For ramjet applications, the nuclear 
reactor offers the possibility of relatively unlimited range or endur- 
ance at low as well as high altitudes. 

In this thesis, the general equations are derived; and two exam- 
ples, a rocket using hydrogen as the working fluid and a ramjet, are 
given. The results of these examples show a distinct advantage in ori- 
ficing to increase the specific impulse of a rocket, but the advantages 
for a ramjet are so slight that it is doubtful if the additional complexity 


in reactor construction arising from orlficing is warranted. 











a 
ii. ANALYSIS AND GENZTRAL SOLUTION 


- ee ee 

In most nuclear reactors the flow of the coolant through the 
cooling passages is turbulent. For this reason the motion of the fluid 
cannot be defined exactly, 30 exact analytical solution of flow and heat 
transfer problems is impossible. However, using the equations of 
motion, where the laminar transport properties are replaced by ef- 
fective values which account for the effects of turbulence and which 
must be determined experimentally, we are able to solve these prob- 
lems. 

Tne Reynolds analogy, first postulated by Osborne Reynolds in 
1890 and verified by experiment, states that momentum and energy 
are transferred in the game way in turbulent shear flow. This anal- 
ogy can be expressed quantitatively by the equation: 


Ke Dees (2. 1) 





where qo is the heat flow rate per unit area normal to the surface 
and T, is the shear stress at the wall. For fully developed turbu- 
lent flow, the velocity profile normal to the surface is nearly flat; 
hence, integrating between the wall where u is zero and a station 


where u and a can be considered to be mean values 


$e ” ie 


BENT =1,) 0 pu 








(2. 2) 


or 


CG, 7 44 (2.34 








rae 


where Cy and oF. are the heat transfer coefficient (Stanton number) 
and friction coefficient respectively. 

If we further write the equations of momentum and energy, 
averaged with respect to time, for mean flow parallel to a surface we 


find: 


6 


Use -puv' + A/a (2. 4) 


4 


Fo 


(2. 5) 


where B is the viscosity and k is the conductivity, and upon re- 


arrangement to better show the similarity: 


oe ee p) Aw 
t, = p ( Au/dy a Ay 
= , (2. 6) 

aes Oy + ee) a 
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where ) = fp is the kinematic viscosity and  ¢ = fa is the Prandtl 
number. Thus, it can be seen that if the Reynolds analogy holds for 
the turbulent properties, and if 7~= 1, Equation 2.3 holds through- 
out turbulent and laminar regions. 

For most gases g¢ —» 0.7 - 1.90, so the Reynolds analogy can 
be assumed to hold reasonably well. The analysis and examples be- 


low are restricted to gas~-cooled reactors with turbulent flow in the 


coolant channels. 











é.2 Temperature Distribution 

Due to non-uniformity of the heat generation in the core of the 
reactor, the average temperature rise of the coolant through the 
reactor core must be determined to find the overall enthalpy rise 
through the reactor. The temperature rise in the individual channels 
must be found first. Integration of each coolant channel temperature 
rise over the reactor cross-section gives the total enthalpy rise. 

The heat balance for a length of coolant passage, dx = /dé, at 


a given radial station, fi » is given by: 
. T= ge AL f(s) aa (2.7) 
Mm Cp Adal , F(S 4+ © 


where m is the mass flow rate per unit area, oe is the specific 

heat at constant pressure, A is the cross-sectional area of the 
coolant passage, T is the static temperature of the coolant, qo is 
the power density at the point of maximum heat generation in the 
reactor core (usually at the center of the core), & = 7 is the dimen- 
sionless coolant passage length measured from the channel entrance, 

» = is the dimensionless radius of the core measured from the cen- 
ter, {(€) is the power distribution function in the axial direction, g(y ) 
is the function in the radial direction, and « is the porosity, or void 
fraction, defined as the ratio of the total coolant passage cross-section- 
al area to the total cross-sectional area. In general, “ could bea 
function of » . 


Rewrite Equation 2.7 as: 


a - AT, t(§) (2. 7a) 


aN 








2G. 


oe ( 

where \ = aie dal is a dimensionless parameter. Since m 
Puce Ts of 

is constant for a given coolant passage, Equation 2. 7a can be integrated 


Lo. 


ae Te + din (fierasl (2. 8) 


The heat balance between the wall and the coolant is given by 


the relation: 


mcpAaT = Crheuce(Ty-T) Pldg (2.9) 


where CL is the Stanton number anc depends chiefly upon the Rey- 
nolds number for gases with Prandtl number close to unity; p is the 
coolant density, u is the velocity of the coolant, T is the channel 
surface temperature, and F is the wetted perimeter. 

An empirical relation used in flow through pipes of circular 
cross-section is: 


_ O.O2.5 
C, = 


}. oO. 2.190 
Re aa G { ) 


where Re is the Reynolds number. The Reynolds number is de- 


fined as: 


a ji eS (2.40 


where 4A/P =D, , the hydraulic diameter. Since R, for a constant 


h 9 
area channel varies only with be which in turn depends cnly on the 
temperature and is approximately proportional to the three-fourths 
power of the absolute temperature, the variation of Gy (and of C €! 


the friction coefficient) along the passage will be quite small, partic- 


ularly for turbulent flow. 
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Since m= pu, Equation 2. 9 can be written as: 


AT _ 7 
j. a(t, aT) (2. 12) 


where 6 = Cy a , ig constant, and is a dimensionless parameter 
proportional to the heat transfer coefficient, or Hence, the wall 


temperature is given by: 


Teme ay Ge (2. 13) 
or, from Equation 2. 8: 
5 
ee T + AT ASE) + fF (gidst (2. 14) 
aYw _ 


The masdmum wall temperature is given by the condition, aa ae 


> 


A =o 
ae | $45) + at 3) (2. 14a) 


From this equation, the vaiue of € for the maximum wall temperature 
can be found, and substituting into Equation 2. 14 TY, » the maximum 
wall temperature,can be obtained. . 
2.3 Pressure Distribution 

The coolant will undergo a pressure drop as it flows through 
the channels. This pressure drop arises from two factors; the ac- 
celeration of the coolant as it passes through the channels, and the 
friction drag opposing the flow of the coolant through the passages. 
The acceleration cf the flow results from decreasing density of the 
coolant, which in turn results from increasing temperature and de- 
creasing pressure. The friction drag cannot be avoided since the 


Reynolds analogy holds. 





ae 


The pressure drop along an individual channel can be determined 
by applying the equations of momentum and continuity. The momentum 


balance for an elernent of length of the coolant passage is given by: 


8 eu. wei+tGpu Bsc (2. 15) 


AlF 


where p is the static pressure of the coolant and Ce i6 an empiri-~ 
cal friction coefficient. From the continuity and perfect gas relations: 


i ngGahame™ T.(; t)=¢ ela ae ~ pete 1) 


where R is the gas constant and subscripts, o, refer to the entrance 


BP 16 


Station of the coolant channel. Substituting Equation into Equa- 


tion * 15 and eu 


=. 2 
1) dp us (ong T| =o 2. 1% 
[p-Serclae+ der (b+ a 
where Be = = é: , a dimensionless parameter proportional to Ce 


the friction coefficient. 

The second term in the coefficient of dp/d& in Equation 2 17 
is 2 result of the acceleration due to the pressure drop in the channel. 
In many applications this term can be neglected and then Equation 2.1% 
can be integrated exactly; however, in flow through a reactor for pro- 
pulsion applications, the pressure drop may be quite large, so this 
term must be retained. The term containing dT/d& is a result of 
acccleration due to the heat input to the coolant as it flows through 
the passage, and will be quite large for nuclear reactors. Equation 


é. 17 is an exact equation for the pressure drop, but it cannot be 


integrated in closed form. However, the equation can be integrated 








alee 
by an iteration process if the Mach number is low. 
2.3.1 Iteration Procedure 
Introduce a dimensionless parameter, € , defined as: 


€ = (rn) Re = (P.ue) RT a os 


< 
Pp. ent ar «! 





The subscripts, o, refer to the channel entrance, and since the 
product RT, igs proportional to the square of the local speed of 
sound for perfect gases, € is seen to be ye’, and thus propor- 
tional to the a number at the channel entrance. It will be as- 


sumed that 7S ee) a te <«<!| as well. 


Equation 2.17 becomes, upon substitution of € , 


£ -¢\ 4p fr egal are a (2. 18) 
Pe P/. ite a T, de ey} 
Let p be expanded in a power series in é€ : 
of £ Cz 2, 
Po opast = ee eee Sa) 
Substituting this series into ae 2.18, we optain: 
fr) (2) 
[ + a wy 7 Jam |] Ae wie (£-| re | (2. 20) 
Te ae AE \ Pe 


a : 
ae 


and, upon collecting powers of € , the differential equation for pi}? 


is found as: 


: 2. 21 
2. (=) +t ( 30 + &T) ve (2. 21) 





21 Sia 
and integrating Equation 2. 2l, 
fv) aes 


lon 
es, Te + & T) 48 —_ 


» 
For the second approximation, the equation for pi4} is: 


a f ,f) es tJ 
ag (i]t (R-%) (Eo) =e a 


and integrating Equation 2. <3: 





5) 3 a) . (toy 
ico. 2 ((e-0) 4 £-) as 


2 a. ok 
P. Po Ede | Be ( ) 


Define the static pressure loss as: 


ae be 7 (p+ Be) Be eye tees. ) 


(> 
where is the first approximation for p at station 1, the reac- 
¢ 
tor exit, and P is the second approximation at this station. There- 
fore, the pressure loss ratio may be expressed ag: 


ees ie 


2 = re 
i 
= pen ((4t - AT) As (2. 25) 
rl Soe ' , 
t e*{—t, ar oAT{ (UG wpe tugs THY AF 
a ‘ 


© wae 


or 


(2. 26) 


These functions, Fy and Fo, may be expressed in a inore con- 


venient form by integration and rearrangement to the iollowing: 





=i 42 


- = 3, Al Fas + a {8 #1545 | 
| . 


(2. 27) 


ee a 
eee (qe 4 (7 . (sgp45| 


‘a - ee -8 

+ < — 8 f ? 

wate BT (T-TJdE the OTl era eat yet a (2. 28) 
aes a: Te {3 oe ote 2 5 


2.4 First Crder Approximation 
The function, f{€), has been defined as the axial power density 


distribution. Hence, Fy may be written as: 


ee +A[ B+ 4 (8-8) (2. 29) 

where Bo and By are defined as: 
ios S £(§)45 = const. (2. 30) 
B= [ EF(§)Agy = const (2. 31) 


If we further require that £(€) be symmetrical about the center of 
the reactor, § = 3, integration of B, shows that B, = 3B,: There- 


fore, Fy », under these assumptions becomes: 


Fo= (3¢ tA BL (it MY ) (2. 32) 


In order to integrate Equation 2. 22 to find the first order ap- 


proximation to the static pressure loss across the reactor, Equation 








a6 


2.7 can be rewritten as: 


Se AT F(S)= KT aly — Fi8) (2. 33) 


where K is defined as: K = ee! 


and the parameter, A , becomes: 


: (mH) 
A: KA —— 


The exit temperature of the working fluid as a function of radial 
position can be determined by integrating Equation 2.33, since the 


integration is in the axial direction only: 


=e ee 2) 2. 34 
T() -T KT, 2 a (2. 34) 
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Using these results, Equation 2. 26 for the first approximation to the 


pressure loss can be expressed in either of the following equations: 


Ap ie pe 
ais € 3, ™ KBB (1+ By) Ve + (2. 35) 


or, alternatively, using the relation of Equation 2. 34: 
ee a (3¢ + Ty = ( + Bef) € amet (2. 36) 


Neither of the two Equations 2.35 and 2. 36 is in a convenient form 
for integration. However, without changing the value of the equation, 


— To 
we may add and subtract the term, € a (+ (34/,. ) , trorn the 
a, =—To 


right side of the equation. is at present unknown and is 


tc be determined. Ferforming this operation, Equation 2.35 becomes: 


OP = = e(pe+ THE fe ef) | e(i+ Pes )(KB$ 4vé-L é| (2. 37) 





*l6= 


—p 


For a given reactor under given flow conditions, Jigen will be a 


— 


constant. Further, when integrated over the radial area of the reac- 
tor, T, and T, will not differ to a great degree, so the second 
term on the right side of Equation 2. 37 will now be very small when 
compared to either of the other two terms. 


For convenience, rewrite Equation 2. 37 as: 


Ap . cc, vy ok Rees 
a ~ Five) + | cve (ve | (2. 38) 


or, rearranging in terms of vo S 


J]. dpe 
(Ve) + [eve - Bey 7 Fees aed (2. 39) 


where: Ve 


(2. 49) 





Since the middle term in Equation 2.39 has been purposely made 
very small in comparison with the other two terms, this equation can 
now be solved by successive approximations. Write Equaticn 2.39 as: 
Lim *, 5(£e - Bey" - ore} - 
i (1é) a 5( & Ve E ( ‘all = 


(2. 41) 


and assume a golution of the form: 


fe : -. a,( + b, 3 +b, 3 Ae )| (2. 42) 
> 


After substituting this assumed solution into Equation 2. 41, col- 


lecting similar powers of oS , and taking the limit as a> |, ae 





oie 


obtain: 


———= 
_ APP. | | G D.| 
Conse ip ee Lagos 2 is 
aad - . ° Ex E 


e 2. 
++ = - SER + agin 
8 = 4PM. EVE Sr, E 


(2. 43) 





It is more compact to leave Equation 2. 43 in terms of C, D, and E, 
since they are defined in Equation 2.40. For integration in the radial 
direction, it can be seen that C contains the variable of integration, 
g( Y) ). For the first approximation, use the first term (which does 
not contain C or D); similarly, for the second approximation to ie 
terms linear in C and D only are required. Higher order requires 


higher powers, 


2.5 Average Temperature Rise Through Reactor 


The total enthalpy rate of increase, AQ, in the exhausting gases 


is: 


< p(t, -T.) amram {a mr 9 ay (2. 44) 


| 


\ 
co (T-L)arep/ee ic Je iat 

where m is defined as the total mags flow rate found by integrating 
the individual coolant mass flow rates over the cross-section of the 
reactor, the integrated average value of os is used, and a is the 
radius of the reactor. 


However, AQ may be found directly as: 





~ ps. 
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21a cp j m(n; Toy “Td 


AQ - 
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4 
samara CB | gin ie os 


Therefore, by equating Equations 2. 44 and 2. 45 we may solve for the 


average temperature rise through the reactor 


= 7 . hs: Sg) yay 


cp PVRT. fxim) Ve(m 44 
| (2. 46) 
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The exit temperature at the centerline of the reactor is found 
by putting f equal to zero in Equation 2, 34 


T (o\- ly 


il 
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em 
FES 


(2. 47) 
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ie 4, be 


(2. 478) 
Hence, dividing Equation 2. 47 by 2. 46, the ratio between the center- 


line exit temperature of the reactor and the average esdt temperature 
is: 


Tel -T. gro) ee 
——— 


(2. 48) 
aoe “ee Elo) {gin yan 





It is not convenient to express the relation between the maximum 


wall temperature and the average temperature increase across the 
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reactor in general. However, if f(§) amd g(7) are known functions, 
this relation can be found quite readily as will be shown in the following 


examples. 
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HI. EXAMPLES 


3,1 Reactor Calculation Using Specific Power Distribution Functions 
Thus far in this thesis the general equations for determining the 
overall temperature rise and pressure drop through the reactor have 
been obtained while leaving the specific power distribution as general 
functions. An example wiil now be given showing how the results ob- 
tained in the general derivation might be applied to a given reactor. 
The general solutions for the neutron density distribution in a 
cylindrical nuclear reactor are given in terms of trigonometric and 
Bessel functions. However, for convenience, a parabolic neutron 
density was assumed since the integrated area under a parabola dif- 
fers but slightly from the more generally accepted solutions of Bessel 
functions or trigonometric functions. For this example, assume a 
rignt-circular cylindrical reactor having parabolic normalized power 


distribution characteristics as: 


FCs) = | =A (3.1) 
qq) = | (3. 2) 


Further, assume the porosity ie uniform across the reactor. The 
power distribution paramcters to be perturbated will be @ and 4, 
which determine the relative ''flatness"’ of the asdal and radial power 
input cf the reactor. These parameters are primarily a function of 
the amount of reflector material used in the reactor design if homo- 


geneous fuel loading is assumed. 
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With these functions, Bo is found from Equation 2. 30 to be: 


{F(sas ae eZ (3. 3) 


Using Equation 2. 48 and working out the integrals indicated to 


the second appro:dmation in /é€ and with: . 
glo) = zs - 
— 4) )A + Salt eZ) Cm | 4) 
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we find, after simplification, that: 
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oa = = aa = 3. 5) 
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To find the maximum wall temperature, which will occur on the 
| centerline of the reactor, use Equation 2. 14a to determine the value 
of & to put in Equation 2.14. Fram Equation 2. l4a, tT, occurs 
when 

= eae mes + J 

S = 77 6 Yah" 46 


and hence, from Ecuation 2. 14: 


i t. 1) eet = 2 
< + Nf ($8 342 24+ a AO 7 2 | (3. 6) 





and using Equation 2, 47a: 
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Hence, given a maximum wall temperature limitation in the 
reactor, the relation (as a function of 9, 6, £, and Bs) between the 
average temperature increase of the gas flowing through the reactor 


may be found by multiplying Equations 3.7 and 3.5: 
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It is simpler to leave Equation 3.8 in this form; hawever, for given 
values of the parameters of the e ‘vation, it reduces to a simple 


yeahs for a given Swern a To 
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quadratic to determine 
If all wetted surfaces of the channel are heat transfer ©. exteCee 
and from the Reynolds analogy, it is reasonable to assume that f = 8 f 
for fully developed turbulent flow. Furthermore, it is assumed that 
= § for comparison purposes, though neither of these assumptions 
are requirements of the analysis. The results obtained from using 
Equation 3.8 are shown in Table l. This table shows the maximum 
wall temperature required to obtain a temperature increase of 
1000°R and 1500°R for an initial gas temperature of 1000°R for var- 
ious values of Band $6. The table was prepared on the assumption 
that ¢ = 6 andf = Be . The required temperature increases were 
selected quite arbitrarily, but they furnish an indication of the trend 
in masximurn wall temperature versus average temperature increase 
through the reactor. From Table I, it can be seen that the maximum 
wali ternperature required for a given temperature increase of the 


working gas is markedly increased as the power distribution va ies 











meas 
away from a linear relation (46 equals zero) or when the coefficient 
of heat, ff , gets very low. 


3.1.1 No-orifice Condition 
The static pressure loss may be found by using Equations 2. 44 


aad 2. 45: 





Taken AR 
= 2 Tae. (Ri KT Bf giqinan 
(3. 9) 
Or Mm T= tk i eo VRT Cee ca 
fea Ue 


But, from the simplification procedure (not shown) used in deter- 


mining Equation 3.5, it was found that 
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therefore, Equation 3.9 can be expressed as: 
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towever, 1 
xo Ta 
Bo Eiquation 3. 10 can be reduced further to: 


is just the mass flow rate, PU, Per channel, 
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Figure 1 shows the variation in static pressure loss with 


VE (which is proportional to the mass flow) for various values of 


Be for an assumed temperature rise ratio of unity across the reactor 
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3.1.2 Artificial Flow Restriction (Orifices) 

For installed orifices scaled so that the mass flow in the outer 
channels is restricted to the correct value to give a uniform tempera- 
ture rise across the reactor, the average exit temperature will be 
identical with the exit temperature from the center channel of the un- 


restricted reactor. Therefore, from Equation 3.7 








tay te ( Ta ; ( ie: 
: 3 Se a -¢| (3. 12) 


If the orifices restrict the flow correctly, the pressure decre- 
ment they provide will make all channel pressure losses identical 
with the pressure loss in the centerline channel, which {5 unrestrict- 
ed. Therefore, the first approximation to the pressure loss ratio 
can be determined from Fquation 3. 4 by assigning a value to /E(o) 
under given reactor conditions. 

The other change that will result from artificially restricting 
the channel entrances will be the total mass flow. This total mass 
flow can be determined directly using Equation 3.11. Since the 
average temperature of the gases is the exit temperature, a constant 
across the reactor, by using Equation 2.34 we are able to determine 
m if it is recalled that it is a factor in the parameter, Jé . Hence, 


integrating, the total mass flow is: | 


m= 274 [am dy ae } fT 8 


ee 
As 


eee 


a poeny lialg 
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However, since |, re, is the exit temperature of the centerline chan- 


nel, Equation 3.12 can be rewritten by using Equation 2. 47 and the 
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3.2 Rocket Example 

Using the assumed reactor power distribution characteristics 
of Equations 3.1 and 3.2, a rocket engine using molecular hydrogen 
aus the vorking fluid was analyzed. No dissociation effects were con- 
sidered. The fuel was assumed to be liquified in the propellant tanks; 
however, due to cocling requirements of the nozzle and reactor reflec- 
tor, the hydrogen was assumed to have reached a gaseous state before 
entering the pumps of the "combustion chamber". 

Account was taken of the power required to drive the pumps by 
using = portion of the gas from the reactor exit to drive a turbine 
which furnished power for the pumps. The temperature of the gases 
as they were bled frorn the reactor exit was reduced to acceptable 
turbine operating temperatures by mixing it with the required amount 
of coolant hydrogen to lower the turbine gas inlet temperature to an 
assumed value of 2000°R. oe 
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For this example, define the following parameters: 
Mp mass flow rate of pump 
Mp mass flow rate of nozzle 


3 bleed fraction to drive auxiliary turbine, therefore, 
Mp = Me C1 ts } ) 


‘Ig reactor exit-entry total pressure ratio = pi, /5/_ 


Ty assumed exit-entry total pressure ratio due to mixing 
loss from exit temperature gradient = p3_ 7 


T, nozzle exit-entry total pressure ratio (assumed }. 0) 


The thrust, F, of a rocket is given by the equation: 


EF o= Me Us + Ag (pa - Pa) (3. 14) 


where Pr is the ambient pressure and subscript 4 refers to exit 


station. The specific impulse, I, is found by the equation: 


Bs ollie Pg - Pa) 
m,4 g ae Pa (3. 15) 


Using the defined parameters, the nozzle mass flow rate, Me ‘ 
can be found using the continuity equation and the perfect gas relations 


for a choked throat as: 








e +! 
mM ae | Z (5-1) 
Ne Pry = ml Vie e+ (3. 16) 
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lf the expansion ratio, Pat ‘p ,, 18 specified, the nozzle-exit 
ratio can be determined from the following equation if isentropic 


flow ia assumed in the nozzle: 
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Me (PA) TL ~ (rg) ‘ (3. 17) 





where rary = itm P ie . Equation 3.17 uniquely defines the 


area ratio and subsequently the exit Mach number from the relation: 


St! 

2( 8-0) 
a 4 pec M> ) (3. 18) 
Ay M, Sti e 


Hence, M i is now known. 





Since M 4 has now been determined, the exit velocity, u 4 can 


be found if T, is known. 
— 


U, ~ Ta, (1+ sy) 
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Therefore; ae ' : 3.1 

| TU (tts (re ge) a 
where a is the static ternperature increase found from Equa- 
Peon 6.5, T) is assumed equal to Tyo since the entering Mach num- 


5 
ber will be very small, and M, was determined from Fig. 6, Dailey 


and Wood (2), modified to account for the friction in the coolant pas- 
sages. 


Knowing Ty T4 is now known, and the exit velocity, Uys is: 
oa 
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The total pressure ratio across the reactor is found in a similar 


manner using the perfect gas relations: ‘ 


Ap at 2\ «6-1 
i ae eee — = (1 = Sp.)/i4 my) 
R ee Fa 2 





(3. 22) 


where - is the static pressure ratio across the reactor as deter- 
( 
mined from Equation 3. 10a by assigning various values of €. After 


having determined Te? the static pressure of the exhaust is given by 


_ 


the relation: eon 
fal =| a- | 
a ree oa 


ail eae 3: 22 
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The total mass flow rate that must be pumped by the pumps is 


known if € is assigned, since from Equations 3. 10 and 3. ila: 


ee Bx ot 
Tat P. (eT (3, 23) 





In order to pump the required total amount of mass flow, the pumps 
expend power which was assumed to be supplied by the turbine. Equa- 


ting the turbine power output to the required pumping power: 


ees cp A Teure, Mrurb, 
e \p We (3. 24) 


Therefore, assuming that the mixture of bled gas and coolant gas 
enters the turbine at 2000°R and exits at the reactor entrance tem- 


perature, the amount of bled gas from the reactor exit is given by: 


. 5 Luma ic. 
aM, = Pcp p (3. 25) 





-~ 28% 


In working out the example, a "rubber'' engine was assumed, 
i.c., the rocket was assumed to be correctly designed for each vari- 
ation in the parameters. A graphitic reactor is used as the basis for 
assuming the high maximum wall temperature (Ref. 3) since it was 
visualized that a rocket engine would not have a long operating life. 
No account was taken of any possible mass addition to the flow from 
the reactor material itself. 


Assumed conditions 


P= 50 atmospheres 

= 0.1 atmoapheres ( 55, 000 feet) 
° 

Te.) 5500 R 

BP 0.90 for no-orifice condition 


1.00 for orifice condition 


XO =a 0.5 

T)= 500°R 

[21/p,= 50 

yo = yes 

C, = 3,8 Btu/Ib °R 


Using these conditions, the exit-throat areca ratio is 5.55 and the 
exit Mach number is 3.17. 

Figure 2 shows the specific impulse versus 6 obtained from the 
calculations for ¢= 9 = 0.5. A linear power distribution (6 = zero) 
is also shown for comparison. No attempt is made to distinguish 
between the values of mass flow used at each 6, since their effect on 
the specific impulses was negligible, although their values do affect 
the pressure loss across the reactor. From these results, the spe-~- 


cific impulse appears to be primarily a function of the temperature 
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rise imparted to the gas in passing through the reactor. Therefore, 
to get high specific impulses in a rocket, it would seem preferable to 
attempt to increase the temperature of the coolant gas by increasing 
the heat transfer coefficient, although this lowers the allowable mass 
flow rate to preclude very high pressure losses, or choked flow, in 
the coolant passages due to the higher heat input to the gases. Yor 
this same reason, orificing appears tc be advantageous in increasing 


the specific impulses. 


aS a ee 

In the ramjet example, the conventional burner is replaced by a 
nuclear reactor having the power distribution characteristics of Equa- 
tions 3.1 and3.2. This analysis was based on standard ICAO condi- 
tions at 1000 feet altitude. 

In working out the example, a "rubber'' engine was assumed and 
the nozzle was agsumed to be designed to expand to atmospheric pres- 
sure unless tne throat static pressure was lower than ambient. In 
such a case, no expansion was considered, and the nozzle was assumed 


to be terminated at the throat. 
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For this analysis, define the following parameters; 
\T, vatio of free-stream total preasure to static = Por/p, 
iT,  vatio of total pressure across diffuser = A 
Te ratio of total pressure across Reactor = Psr/p,_ 
Tm ratio of total pressures due to assumed mixing loss 
in non-restricted case = Bon p ee {Assumed | for 


restricted case. } 


TR, ratic of nozzle total pressures (assumed 1) 


The thrust, F, of a ramjet is given by: 


ce rau, (4s ~1) t Asps 1 - Pes.) (3. 26) 


Yue to the fact that ramjet engines must operate for long peri- 
ods of tirme, and in an oxidizing atmosphere, a diiferent type reactor 
material must be used in a ramjet thas in a rocket. Because of this, 
the materials suitable for ramjet applications can withstand a much 
lower maximum wall temperature than those for a rocket engine. The 
heat input to the coolant gases will be much less, and consequently the 
acceleration of the gages in the coolant passages will be less. There- 
fore, in this example, it was assumed that the static temperature 
rise and the static pressure loss across the reactor were approximately 
equal to their total, respectively. 

4s4, can be found by finding the ratios of exit Mach number to 
the free stream Mach number and exit total temperature to the initial 


total temperature. 
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and from this we obtain: 
Zz oY z 
ocean q-4 MA 
M- = ES TT, Te ™, : | : = = ae | (3 29) 
M, Ps see! Age ra MA 
2 “ 2 . 
Aliso: 
Ws) 2 Use 2), AT 
To, nee ee (3.30) 
AT : tr -T. ; 
where Te ee of the general reactor analysis. 
Therefore, the exit velocity can be expressed ag 
ie a | _ _ ae) a 
( Or fms or Tents ae (3. 32) 


Since it was assumed that the static pressure loss and total 


pressure loss are approximately equal, Tr is found by using Equa- 


tion 3. 10a to determine Ap/p for various values of €. 


Ee ies 
ta oa (3. 32) 
Knoving “pe the throat pressure can be found from: 
os 
a z o-{ 
Pope = (Ger mR Me Oe 


If this ratio is less than unity, no further expansion was assumed; 


if the ratio was greater than unity, ideal expansion to ambient pressure 


was agsumed. 
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Figures 3 and 4 show the results obtained for Ce (based on the 
reactor cross-sectional area), the thrust coefficient, versus Mach 
number for various values of 6 at an assumed maximum reactor wall 


temperature af 2960°R. C_. is defined as the total thrust divided by 


EF 


the free stream dynamic pressure and the cross-sectional area of the 


veactor. Values of a used in the example were: 


M cP 
1.5 —-25 0.85 
DS oa ES oO 7(S 
Ses 2.6 


ws the assumed mixing loss, was arbitrarily assigned a value of 
0. 90 for the non-vestricted example and unity for the orificed exam- 
ple. Results shown are for an assumed reactor porosity of 0. 5. 

The curves shown are for € = 0,03, but other massa flow 
rates give similar curves. The one shown is about optimum, since 
it was found that further increase in mass flow rate increased the 
pressure loss across the reactor to a point where the overall thrust 
is reduced. While it should be stressed that the initial conditions and 
reactor configuration chosen were arbitrary and this example was 
worked to demonstrate the methods, still, the trend of the results 
should be realistic, if not their absolute magnitudes. From these 
two figures at the arbitrarily chosen configuration, optimum per- 
formance is obtained at 68 of approximately 3-5 at a Mach number of 
about 3. Figure 4 shows that while Ce at the optimum Mach number 
is slightly increased by orificing, the vehicle performance drops 
more sharply away from the optimum Mach number, go it is doubtful 


that there is much advantage to orificing for a ramjet. 
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TASLE f 


Masdmum Wall Temperature Required to Obtain Given Average 
Temperature Increase of 


Coolant Gas (Initial gas temperature assumed to be 1000°R) 


oo 











AT = 1000 = AT = 1500 

p = Pr T yf) —_ 
m 
20 3430 3440 
10 3460 3470 

b= 1.0 

5 3520 3625 
1 3860 7000 
20 2790 2818 
10 2801 2868 
5 2828 2981 
i 2950 4760 
20 2418 2518 
10 2420 2533 
#2 .50 5 2430 2620 
l 2480 3720 
20 2172 2450 
10 2175 2452 
5 2177 2490 
l 2193 3190 
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TABLE I 
(continued) 
AT = 19005 AT = 15900 
Be Be Te) Tw Ti) Ty 
mM 
20 2000 | 2000 2500 2590 
10 2000 2900 2500 2509 
5 2000 2000 2500 2500 


) 2000 2006 2500 2500 
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FIG.@—- SPECIFIC IMPULSE VS. HEAT COEFFICIENT FOR 
NUCLEAR- HYDROGEN ROCKET EXAMPLE 
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FIG.3 — RAMJET EXAMPLE (UNRESTRICTED COOLANT 
PASSAGES) — THRUST COEFFICIENT VSapRee 
STREAM MACH NUMBER 
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FIG.4 - RAMJET EXAMPLE (ORIFICED GOOLANT 
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